home *** CD-ROM | disk | FTP | other *** search
/ AmigActive 26 / AACD 26.iso / AACD / Programming / ace_gpl_release / src / dp / dp_main.s < prev    next >
Encoding:
Text File  |  1996-08-28  |  11.2 KB  |  693 lines
  1. ;
  2. ; Double-precision math interface library for ACE.
  3. ; Copyright (C) 1998 David Benn
  4. ; 
  5. ; This program is free software; you can redistribute it and/or
  6. ; modify it under the terms of the GNU General Public License
  7. ; as published by the Free Software Foundation; either version 2
  8. ; of the License, or (at your option) any later version.
  9. ;
  10. ; This program is distributed in the hope that it will be useful,
  11. ; but WITHOUT ANY WARRANTY; without even the implied warranty of
  12. ; MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  13. ; GNU General Public License for more details.
  14. ;
  15. ; You should have received a copy of the GNU General Public License
  16. ; along with this program; if not, write to the Free Software
  17. ; Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.
  18. ;
  19. ; Register a6 is preserved and restored by each routine
  20. ; since Sozobon C uses this as a link register.
  21. ; 
  22. ; Author: David J Benn
  23. ;   Date: 11th February 1996,
  24. ;      16th August 1996,
  25. ;      29th September 1996,
  26. ;      1st,18th October 1996
  27. ;
  28.  
  29. TRUE    equ    -1
  30. FALSE    equ    0
  31.  
  32.     ; Basic functions.
  33.     xref    _LVOIEEEDPAbs
  34.     xref    _LVOIEEEDPAdd
  35.     xref    _LVOIEEEDPCeil
  36.     xref    _LVOIEEEDPCmp
  37.     xref    _LVOIEEEDPDiv
  38.     xref    _LVOIEEEDPFix
  39.     xref    _LVOIEEEDPFloor
  40.     xref    _LVOIEEEDPFlt
  41.     xref    _LVOIEEEDPMul
  42.     xref    _LVOIEEEDPNeg
  43.     xref    _LVOIEEEDPSub
  44.     xref    _LVOIEEEDPTst
  45.  
  46.     ; Transcendental functions.
  47.     xref    _LVOIEEEDPAsin
  48.     xref    _LVOIEEEDPAcos
  49.     xref    _LVOIEEEDPAtan
  50.     xref    _LVOIEEEDPSin
  51.     xref    _LVOIEEEDPCos
  52.     xref    _LVOIEEEDPTan
  53.     xref    _LVOIEEEDPSincos
  54.     xref    _LVOIEEEDPSinh
  55.     xref    _LVOIEEEDPCosh
  56.     xref    _LVOIEEEDPTanh
  57.     xref    _LVOIEEEDPExp
  58.     xref    _LVOIEEEDPFieee
  59.     xref    _LVOIEEEDPLog
  60.     xref    _LVOIEEEDPLog10
  61.     xref    _LVOIEEEDPPow
  62.     xref    _LVOIEEEDPSqrt
  63.     xref    _LVOIEEEDPTieee
  64.  
  65.     ; Miscellaneous external references.
  66.     xref    _AbsExecBase
  67.     xref    _LVOOpenLibrary
  68.     xref    _LVOCloseLibrary
  69.  
  70.     xref    _MathTransBase
  71.     xref    _LVOSPFieee    
  72.     xref    _LVOSPTieee    
  73.  
  74.     ; Exports.
  75.     xdef    _MathIeeeDoubBasBase
  76.     xdef    _MathIeeeDoubTransBase
  77.  
  78.     section    code, dp_code
  79.  
  80. ;************************************
  81.  
  82.     PUBLIC    _dp_open
  83.     ;
  84.     ; Open IEEE DP libraries x 2.
  85.     ;
  86. _dp_open:
  87.     move.l    a6,_a6
  88.  
  89.     movea.l _AbsExecBase,a6        
  90.  
  91.     ; Open basic library.
  92.     movea.l    #_ieeedoubbaslibname,a1
  93.     move.l    #34,d0
  94.     jsr     _LVOOpenLibrary(a6)
  95.     tst.l    d0
  96.     beq.s    _abort_dp_open
  97.     move.l    d0,_MathIeeeDoubBasBase
  98.  
  99.     ; Open transcendental library.
  100.     movea.l    #_ieeedoubtranslibname,a1
  101.     move.l    #34,d0
  102.     jsr     _LVOOpenLibrary(a6)
  103.     tst.l    d0
  104.     beq.s    _abort_dp_open
  105.     move.l    d0,_MathIeeeDoubTransBase
  106.     
  107.     ; Exit with success code (TRUE).
  108.     move.l    #TRUE,d0    ; okay.
  109.     bra.s    _exit_dp_open
  110.  
  111. _abort_dp_open:
  112.     move.l    #FALSE,d0    ; error.
  113.  
  114. _exit_dp_open:
  115.     movea.l    _a6,a6
  116.     
  117.     rts
  118.  
  119. ;************************************
  120.  
  121.     PUBLIC    _dp_close
  122.     ;
  123.     ; Close IEEE DP libraries x 2.
  124.     ;
  125. _dp_close:
  126.     move.l    a6,_a6
  127.  
  128.     tst.l    _MathIeeeDoubBasBase
  129.     beq.s    _exit_dp_close
  130.  
  131.     movea.l _AbsExecBase,a6
  132.  
  133.     move.l  _MathIeeeDoubBasBase,a1
  134.     jsr     _LVOCloseLibrary(a6)
  135.  
  136.     tst.l    _MathIeeeDoubTransBase
  137.     beq.s    _exit_dp_close
  138.  
  139.     move.l  _MathIeeeDoubTransBase,a1
  140.     jsr     _LVOCloseLibrary(a6)
  141.  
  142. _exit_dp_close:
  143.     movea.l    _a6,a6
  144.     
  145.     rts
  146.  
  147. ;************************************
  148.  
  149.     PUBLIC    _dp_longint_to_double
  150.     ;
  151.     ; Convert LONGINT to DOUBLE.
  152.     ;
  153. _dp_longint_to_double:
  154.     move.l    a6,_a6
  155.  
  156.     move.l    8(sp),d0        ; integer value.
  157.  
  158.     movea.l    _MathIeeeDoubBasBase,a6
  159.     jsr    _LVOIEEEDPFlt(a6)
  160.  
  161.     move.l    4(sp),a0        ; address of double result.
  162.     move.l    d0,(a0)+
  163.     move.l    d1,(a0)
  164.  
  165.     movea.l    _a6,a6
  166.     
  167.     rts
  168.  
  169. ;************************************
  170.  
  171.     PUBLIC    _dp_double_to_longint
  172.     ;
  173.     ; Convert DOUBLE to LONGINT.
  174.     ;
  175. _dp_double_to_longint:
  176.     move.l    a6,_a6
  177.  
  178.     move.l    8(sp),a1        ; address of double.
  179.     
  180.     move.l    (a1)+,d0
  181.     move.l    (a1),d1
  182.  
  183.     movea.l    _MathIeeeDoubBasBase,a6
  184.     jsr    _LVOIEEEDPFix(a6)
  185.  
  186.     move.l    4(sp),a0        ; address of integer result.
  187.     move.l    d0,(a0)
  188.  
  189.     movea.l    _a6,a6
  190.  
  191.     rts
  192.  
  193. ;************************************
  194.  
  195.     PUBLIC    _dp_single_to_double
  196.     ;
  197.     ; Convert SINGLE to DOUBLE.
  198.     ;
  199. _dp_single_to_double:
  200.     move.l    a6,_a6
  201.  
  202.     move.l    8(sp),d0        ; single value.
  203.  
  204.     movea.l    _MathTransBase,a6
  205.     jsr    _LVOSPTieee(a6)        ; FFP -> IEEE SP
  206.     movea.l    _MathIeeeDoubTransBase,a6    
  207.     jsr    _LVOIEEEDPFieee(a6)    ; IEEE SP -> IEEE DP
  208.  
  209.     move.l    4(sp),a0        ; address of double result.
  210.     move.l    d0,(a0)+
  211.     move.l    d1,(a0)
  212.  
  213.     movea.l    _a6,a6
  214.     
  215.     rts
  216.  
  217. ;************************************
  218.  
  219.     PUBLIC    _dp_double_to_single
  220.     ;
  221.     ; Convert DOUBLE to SINGLE.
  222.     ;
  223. _dp_double_to_single:
  224.     move.l    a6,_a6
  225.  
  226.     move.l    8(sp),a1        ; address of double.
  227.  
  228.     move.l    (a1)+,d0
  229.     move.l    (a1),d1
  230.  
  231.     movea.l    _MathIeeeDoubTransBase,a6
  232.     jsr    _LVOIEEEDPTieee(a6)    ; IEEE DP -> IEEE SP
  233.     movea.l    _MathTransBase,a6
  234.     jsr    _LVOSPFieee(a6)        ; IEEE SP -> FFP
  235.  
  236.     move.l    4(sp),a0        ; address of single result.
  237.     move.l    d0,(a0)
  238.  
  239.     movea.l    _a6,a6
  240.  
  241.     rts
  242.  
  243. ;************************************
  244.  
  245.     PUBLIC    _dp_add
  246.     ;
  247.     ; Add two DOUBLE values (result = op1 + op2).
  248.     ;
  249. _dp_add:
  250.     move.l    a6,_a6
  251.  
  252.     move.l    12(sp),a2        ; address of operand 2.
  253.     move.l    8(sp),a1        ; address of operand 1.
  254.  
  255.     move.l    (a2)+,d2
  256.     move.l    (a2),d3            
  257.  
  258.     move.l    (a1)+,d0
  259.     move.l    (a1),d1    
  260.  
  261.     movea.l    _MathIeeeDoubBasBase,a6
  262.     jsr    _LVOIEEEDPAdd(a6)
  263.  
  264.     move.l    4(sp),a0        ; address of result.    
  265.     move.l    d0,(a0)+
  266.     move.l    d1,(a0)
  267.  
  268.     movea.l    _a6,a6
  269.  
  270.     rts
  271.  
  272. ;************************************
  273.  
  274.     PUBLIC    _dp_sub
  275.     ;
  276.     ; Subtract two DOUBLE values (result = op1 - op2).
  277.     ;
  278. _dp_sub:
  279.     move.l    a6,_a6
  280.  
  281.     move.l    12(sp),a2        ; address of operand 2.
  282.     move.l    8(sp),a1        ; address of operand 1.
  283.  
  284.     move.l    (a2)+,d2
  285.     move.l    (a2),d3            
  286.  
  287.     move.l    (a1)+,d0
  288.     move.l    (a1),d1    
  289.  
  290.     movea.l    _MathIeeeDoubBasBase,a6
  291.     jsr    _LVOIEEEDPSub(a6)
  292.  
  293.     move.l    4(sp),a0        ; address of result.    
  294.     move.l    d0,(a0)+
  295.     move.l    d1,(a0)
  296.  
  297.     movea.l    _a6,a6
  298.  
  299.     rts
  300.  
  301. ;************************************
  302.  
  303.     PUBLIC    _dp_mul
  304.     ;
  305.     ; Mulitply two DOUBLE values (result = op1 * op2).
  306.     ;
  307. _dp_mul:
  308.     move.l    a6,_a6
  309.  
  310.     move.l    12(sp),a2        ; address of operand 2.
  311.     move.l    8(sp),a1        ; address of operand 1.
  312.  
  313.     move.l    (a2)+,d2
  314.     move.l    (a2),d3            
  315.  
  316.     move.l    (a1)+,d0
  317.     move.l    (a1),d1    
  318.  
  319.     movea.l    _MathIeeeDoubBasBase,a6
  320.     jsr    _LVOIEEEDPMul(a6)
  321.  
  322.     move.l    4(sp),a0        ; address of result.    
  323.     move.l    d0,(a0)+
  324.     move.l    d1,(a0)
  325.  
  326.     movea.l    _a6,a6
  327.  
  328.     rts
  329.  
  330. ;************************************
  331.  
  332.     PUBLIC    _dp_div
  333.     ;
  334.     ; Divide two DOUBLE values (result = op1 / op2).
  335.     ;
  336. _dp_div:
  337.     move.l    a6,_a6
  338.  
  339.     move.l    12(sp),a2        ; address of operand 2.
  340.     move.l    8(sp),a1        ; address of operand 1.
  341.  
  342.     move.l    (a2)+,d2
  343.     move.l    (a2),d3            
  344.  
  345.     move.l    (a1)+,d0
  346.     move.l    (a1),d1    
  347.  
  348.     movea.l    _MathIeeeDoubBasBase,a6
  349.     jsr    _LVOIEEEDPDiv(a6)
  350.  
  351.     move.l    4(sp),a0        ; address of result.    
  352.     move.l    d0,(a0)+
  353.     move.l    d1,(a0)
  354.  
  355.     movea.l    _a6,a6
  356.  
  357.     rts
  358.  
  359. ;************************************
  360.  
  361.     PUBLIC    _dp_cmp
  362.     ;
  363.     ; Compare two DOUBLE values.
  364.     ; Returns an integer:
  365.     ; 
  366.     ;     +1 if op1 > op2
  367.     ;     -1 if op1 < op2
  368.     ;      0 if op1 = op2
  369.     ;
  370. _dp_cmp:
  371.     move.l    a6,_a6
  372.  
  373.     move.l    12(sp),a2        ; address of operand 2.
  374.     move.l    8(sp),a1        ; address of operand 1.
  375.  
  376.     move.l    (a2)+,d2
  377.     move.l    (a2),d3            
  378.  
  379.     move.l    (a1)+,d0
  380.     move.l    (a1),d1    
  381.  
  382.     movea.l    _MathIeeeDoubBasBase,a6
  383.     jsr    _LVOIEEEDPCmp(a6)
  384.  
  385.     move.l    4(sp),a0        ; address of result.    
  386.     move.l    d0,(a0)
  387.  
  388.     movea.l    _a6,a6
  389.  
  390.     rts
  391.  
  392. ;************************************
  393.  
  394.     PUBLIC    _dp_pow
  395.     ;
  396.     ; Raise one DOUBLE to the power of another.
  397.     ; (result = op1 ^ op2).
  398.     ;
  399. _dp_pow:
  400.     move.l    a6,_a6
  401.  
  402.     move.l    12(sp),a2        ; address of operand 2.
  403.     move.l    8(sp),a1        ; address of operand 1.
  404.  
  405.     move.l    (a1)+,d0
  406.     move.l    (a1),d1        
  407.  
  408.     move.l    (a2)+,d2
  409.     move.l    (a2),d3
  410.  
  411.     movea.l    _MathIeeeDoubTransBase,a6
  412.     jsr    _LVOIEEEDPPow(a6)
  413.  
  414.     move.l    4(sp),a0        ; address of result.    
  415.     move.l    d0,(a0)+
  416.     move.l    d1,(a0)
  417.  
  418.     movea.l    _a6,a6
  419.  
  420.     rts
  421.  
  422. ;************************************
  423.  
  424.     PUBLIC    _dp_abs
  425.     ;
  426.     ; Return the absolute value of a DOUBLE.
  427.     ;
  428. _dp_abs:
  429.     move.l    a6,_a6
  430.  
  431.     move.l    8(sp),a1        ; address of operand.
  432.  
  433.     move.l    (a1)+,d0
  434.     move.l    (a1),d1    
  435.  
  436.     movea.l    _MathIeeeDoubBasBase,a6
  437.     jsr    _LVOIEEEDPAbs(a6)
  438.  
  439.     move.l    4(sp),a0        ; address of result.    
  440.     move.l    d0,(a0)+
  441.     move.l    d1,(a0)
  442.  
  443.     movea.l    _a6,a6
  444.  
  445.     rts
  446.  
  447. ;************************************
  448.  
  449.     PUBLIC    _dp_neg
  450.     ;
  451.     ; Return the negated value of a DOUBLE.
  452.     ;
  453. _dp_neg:
  454.     move.l    a6,_a6
  455.  
  456.     move.l    8(sp),a1        ; address of operand.
  457.  
  458.     move.l    (a1)+,d0
  459.     move.l    (a1),d1    
  460.  
  461.     movea.l    _MathIeeeDoubBasBase,a6
  462.     jsr    _LVOIEEEDPNeg(a6)
  463.  
  464.     move.l    4(sp),a0        ; address of result.    
  465.     move.l    d0,(a0)+
  466.     move.l    d1,(a0)
  467.  
  468.     movea.l    _a6,a6
  469.  
  470.     rts
  471.  
  472. ;************************************
  473.  
  474.     PUBLIC    _dp_ceil
  475.     ;
  476.     ; Return the least integer which is >= a DOUBLE.
  477.     ; Result is DOUBLE not LONGINT!
  478.     ;
  479. _dp_ceil:
  480.     move.l    a6,_a6
  481.  
  482.     move.l    8(sp),a1        ; address of operand.
  483.  
  484.     move.l    (a1)+,d0
  485.     move.l    (a1),d1    
  486.  
  487.     movea.l    _MathIeeeDoubBasBase,a6
  488.     jsr    _LVOIEEEDPCeil(a6)
  489.  
  490.     move.l    4(sp),a0        ; address of result.    
  491.     move.l    d0,(a0)+
  492.     move.l    d1,(a0)
  493.  
  494.     movea.l    _a6,a6
  495.  
  496.     rts
  497.  
  498. ;************************************
  499.  
  500.     PUBLIC    _dp_floor
  501.     ;
  502.     ; Return the largest integer which is <= a DOUBLE.
  503.     ; Result is DOUBLE not LONGINT!
  504.     ;
  505. _dp_floor:
  506.     move.l    a6,_a6
  507.  
  508.     move.l    8(sp),a1        ; address of operand.
  509.  
  510.     move.l    (a1)+,d0
  511.     move.l    (a1),d1    
  512.  
  513.     movea.l    _MathIeeeDoubBasBase,a6
  514.     jsr    _LVOIEEEDPFloor(a6)
  515.  
  516.     move.l    4(sp),a0        ; address of result.    
  517.     move.l    d0,(a0)+
  518.     move.l    d1,(a0)
  519.  
  520.     movea.l    _a6,a6
  521.  
  522.     rts
  523.  
  524. ;************************************
  525.  
  526.     PUBLIC    _dp_sin
  527.     ;
  528.     ; Return the sine of a DOUBLE (radians).
  529.     ;
  530. _dp_sin:
  531.     move.l    a6,_a6
  532.  
  533.     move.l    8(sp),a1        ; address of operand.
  534.  
  535.     move.l    (a1)+,d0
  536.     move.l    (a1),d1    
  537.  
  538.     movea.l    _MathIeeeDoubTransBase,a6
  539.     jsr    _LVOIEEEDPSin(a6)
  540.  
  541.     move.l    4(sp),a0        ; address of result.    
  542.     move.l    d0,(a0)+
  543.     move.l    d1,(a0)
  544.  
  545.     movea.l    _a6,a6
  546.  
  547.     rts
  548.  
  549. ;************************************
  550.  
  551.     PUBLIC    _dp_cos
  552.     ;
  553.     ; Return the cosine of a DOUBLE (radians).
  554.     ;
  555. _dp_cos:
  556.     move.l    a6,_a6
  557.  
  558.     move.l    8(sp),a1        ; address of operand.
  559.  
  560.     move.l    (a1)+,d0
  561.     move.l    (a1),d1    
  562.  
  563.     movea.l    _MathIeeeDoubTransBase,a6
  564.     jsr    _LVOIEEEDPCos(a6)
  565.  
  566.     move.l    4(sp),a0        ; address of result.    
  567.     move.l    d0,(a0)+
  568.     move.l    d1,(a0)
  569.  
  570.     movea.l    _a6,a6
  571.  
  572.     rts
  573.  
  574. ;************************************
  575.  
  576.     PUBLIC    _dp_tan
  577.     ;
  578.     ; Return the tangent of a DOUBLE (radians).
  579.     ;
  580. _dp_tan:
  581.     move.l    a6,_a6
  582.  
  583.     move.l    8(sp),a1        ; address of operand.
  584.  
  585.     move.l    (a1)+,d0
  586.     move.l    (a1),d1    
  587.  
  588.     movea.l    _MathIeeeDoubTransBase,a6
  589.     jsr    _LVOIEEEDPTan(a6)
  590.  
  591.     move.l    4(sp),a0        ; address of result.    
  592.     move.l    d0,(a0)+
  593.     move.l    d1,(a0)
  594.  
  595.     movea.l    _a6,a6
  596.  
  597.     rts
  598.  
  599. ;************************************
  600.  
  601.     PUBLIC    _dp_exp
  602.     ;
  603.     ; Return e to the power of a DOUBLE.
  604.     ;
  605. _dp_exp:
  606.     move.l    a6,_a6
  607.  
  608.     move.l    8(sp),a1        ; address of operand.
  609.  
  610.     move.l    (a1)+,d0
  611.     move.l    (a1),d1    
  612.  
  613.     movea.l    _MathIeeeDoubTransBase,a6
  614.     jsr    _LVOIEEEDPExp(a6)
  615.  
  616.     move.l    4(sp),a0        ; address of result.    
  617.     move.l    d0,(a0)+
  618.     move.l    d1,(a0)
  619.  
  620.     movea.l    _a6,a6
  621.  
  622.     rts
  623.  
  624. ;************************************
  625.  
  626.     PUBLIC    _dp_log
  627.     ;
  628.     ; Return the natural logarithm (base e) of a DOUBLE.
  629.     ;
  630. _dp_log:
  631.     move.l    a6,_a6
  632.  
  633.     move.l    8(sp),a1        ; address of operand.
  634.  
  635.     move.l    (a1)+,d0
  636.     move.l    (a1),d1    
  637.  
  638.     movea.l    _MathIeeeDoubTransBase,a6
  639.     jsr    _LVOIEEEDPLog(a6)
  640.  
  641.     move.l    4(sp),a0        ; address of result.    
  642.     move.l    d0,(a0)+
  643.     move.l    d1,(a0)
  644.  
  645.     movea.l    _a6,a6
  646.  
  647.     rts
  648.  
  649. ;************************************
  650.  
  651.     PUBLIC    _dp_sqrt
  652.     ;
  653.     ; Return the square root of a DOUBLE.
  654.     ;
  655. _dp_sqrt:
  656.     move.l    a6,_a6
  657.  
  658.     move.l    8(sp),a1        ; address of operand.
  659.  
  660.     move.l    (a1)+,d0
  661.     move.l    (a1),d1    
  662.  
  663.     movea.l    _MathIeeeDoubTransBase,a6
  664.     jsr    _LVOIEEEDPSqrt(a6)
  665.  
  666.     move.l    4(sp),a0        ; address of result.    
  667.     move.l    d0,(a0)+
  668.     move.l    d1,(a0)
  669.  
  670.     movea.l    _a6,a6
  671.  
  672.     rts
  673.  
  674. ;************************************
  675.         
  676.     section    data, dp_data
  677.  
  678. _ieeedoubbaslibname:    dc.b    'mathieeedoubbas.library',0
  679. _ieeedoubtranslibname:    dc.b    'mathieeedoubtrans.library',0
  680.  
  681. _MathIeeeDoubBasBase:    dc.l     0
  682. _MathIeeeDoubTransBase:    dc.l     0
  683.  
  684. ;************************************
  685.  
  686.     section bss, dp_bss
  687.  
  688. _a6:            ds.l 1
  689.  
  690. ;************************************
  691.  
  692.     END
  693.